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1 Introduction 

Cosmological inflation predicts the initial power spectrum of scalar density fluctuations to 
be close to scale- invariant with just a slight scale dependence. The simplest proposal for 
the shape of the spectrum is to assume a power-law parameterisation in terms of a spec- 
tral amplitude Ag and a spectral index or tilt parameter rig. Although this form has been 
in good agreement with cosmological observations, recent analyses from the Wilkinson Mi- 
crowave Anisotropy Probe satellite (WMAP; [1]), and the Atacama Cosmology Telescope 
(ACT; [2]) have confirmed that the scale invariant (rig = 1) spectrum is now excluded at 3a 
C.L. Similar results are obtained when measurements from the South Pole Telescope (SPT; 
[3]) in combination with the WMAP data are considered. It has also been shown that if a 
running of the spectral index is taken into account, allowing deviations from the power-law 
spectrum, WMAP+ ACT and WMAP+ SPT data show a preference for a negative run- 
ning value at 1.8a C.L. Thus, consideration of models that slightly deviate from the simple 
power-law might be required. There have been several alternatives proposed. Some physi- 
cally motivated models include an exponential large scale cut-off [4], discontinuities in the 
early universe from phase transitions [5], closed universe inflation [6], models in which the 
power spectrum drops to zero below some cut-off scale [7]. Some use observational data to 
constrain an a priori parameterisation or attempt a direct reconstruction using wavelets [8] , 
deconvolution methods [9-11], binning the spectrum into an arbitrary number of band pow- 
ers [12-15], Bayesian reconstruction [16], principle component analysis [17] and minimising 
the level of complexity needed via a cross-validation [18], amongst many others. 

In this paper, we are interested in selecting the preferred shape for the primordial spec- 
trum using the Bayesian evidence as an implementation of Occam's razor: a simpler model 



should be preferred, unless the data require a more sophisticated model. First, we determine 
the structure of the primordial power spectrum using an optimal model-free reconstruction. 
Our approach considers possible deviations from the power-law parameterisation by mod- 
elling the spectrum as a linear interpolation between a set of 'nodes' which can vary in both 
amplitude and ^-position. Within this analysis we have included phenomenological features 
which might be relevant to the description of CMB observations, such as a large scale cut-off, 
a broken spectrum and a spectrum with a possible turn-over at any position in fe-space. The 
reconstruction process is essentially the same binning format used previously by a number of 
authors, however here we allow the data to decide the level of complexity of the model - the 
number of nodes and their optimum position - via the Bayesian evidence. Then, for compar- 
ison, we compute the Bayesian evidence for a set of existing model proposals: a power-law 
parameterisation including both tilt and running parameter, a modified power-law spectrum 
to account for a drop off at large scales and the Lasenby & Doran (LD) model based on 
a closed universe. Finally, for each model we compare its Bayesian evidence and according 
to the Jeffreys guideline we select the best model preferred by current data. In a previous 
paper [19], we have constrained the parameter space which describes the primordial spectrum 
derived from the LD model. In this work, using an optimal model-free reconstruction, the 
shape of the spectrum is determined directly from the data. 

The paper is organised as follows: in the next Section we describe the basic parameter 
estimation and model selection, we also list the datasets and the cosmological parameters 
considered. In Section 3 we carry out the reconstruction for the primordial power spectrum 
and also consider different existing parameterisations suggested to describe the form of the 
spectrum. We show the resulting parameter constraints and the evidence for each worked 
model. Finally, in Section 4 we decide which model provides the best description for current 
observational data and present our conclusions. 

2 Bayesian Inference 

2.1 Parameter Estimation 

Given a model or hypothesis H for some data D, Bayes' theorem tells us how to determine 
the probability distribution of the set of parameters on which the model depends. Bayes' 
theorem states that 

,_, , Pr(D|0,m Fi(@\H) , , 

Pr(eiD.fl)^ 'p;(D|g)' ' ('■» 

where, for brevity, we denote Ft{&\H) = 7r(0) as the prior probability which is upgraded 
through the likelihood Pt(D\&,H) = C{&) when experimental data D are considered. The 
aim for parameter estimation is to obtain the posterior probability Pr(0|D,i?) which rep- 
resents the state of knowledge once we have taken into account the new information. The 
normalisation constant is usually refered to the Bayesian evidence Pr(D\H) = Z. Since this 
quantity is independent of the parameters 0, it is often ignored in parameter estimation but 
plays the central role for model selection. 

Throughout the analysis we compute posterior probabilities for each model in the light 
of temperature and polarisation measurements from the 7-year data release of the Wilkinson 



Microwave Anisotropy Probe (WMAP; [1]). The ACT and SPT data provide similar con- 
straints on the primary cosmological parameters [2, 3], thus since our results are not likely 
to be significantly affected according to which is chosen, and to avoid possible problems 
due to overlap of the sky regions, we have restricted the current analysis to use of just the 
ACT data. In addition to CMB data, we include distance measurements from the Supernova 
Cosmology Project Union 2 compilation (SCP; [20]) and large scale structure data from the 
Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Luminous Red Galaxy (LRG) power 
spectrum [21]. We also impose a Gaussian prior from measurements of the Hubble parameter 
today Hq from the Hubble Space Telescope (HST; [22]) key project. 

We consider purely Gaussian adiabatic scalar perturbations and neglect tensor contri- 
butions. We assume a flat ACDM model ^ specified by the following parameters: the physical 
baryon Ob/i^ and cold dark matter density Qdm^i^ relative to the critical density {h is the 
dimensionless Hubble parameter such that Hq = lOO/i kms~^Mpc~^), 6 is lOOx the ratio of 
the sound horizon to angular diameter distance at last scattering surface and r denotes the 
optical depth at reionisation. Aside from the Sunyaev-Zel'dovich (SZ) amplitude Asz used 
by WMAP analyses, the 148 GHz ACT likelihood incorporates two additional secondary pa- 
rameters: the total Poisson power Ap at I = 3000 and the amplitude of the clustered power 
Ac- The parameters describing the primordial spectra in each model are listed in Section 3, 
together with their flat priors imposed in our Bayesian analysis. 

2.2 Model Selection 

Model selection applies a similar type of analysis to parameter estimation but now at the 
level of models rather than parameters. The key quantity to bear in mind is the Bayesian 
evidence. It balances the complexity of cosmological models and then, naturally, incorporates 
Occam's razor. Let us consider several models or hypotheses Hi, each of them with prior 
probability Pr{Hi). Bayes' theorem for model selection is 



Pr(g,|D) ^ Pr(D|g,) Pr(g,) 

Fv{Hj\T>) FT{'D\Hj)Fr{Hj)' ^ ' ' 

The left-hand side denotes the relative probability of the model given the data, whereas in 
the right-hand side appears the key quantity to compute: the Bayesian evidence for each 
model. This evidence is nothing more than the average of the likelihood over the prior 



f c{e)7T{e)d^@, (2.3) 



where N is the dimensionality of the parameter space. Note that a complex model will usu- 
ally lead to a higher likelihood, but if the fit is nearly as good as a simple model, the evidence 
will favour the simpler model through the smaller prior volume. 

For convenience, the ratio of two evidences Zq/Zi (or equivalently the difference in log 
evidences IhZq — InZi) is often termed the Bayes factor ;Bo,i- 

Bo,i= ln|^. (2.4) 

Zi 

^Except for the LD model, which is based on a closed universe Qk < 0. 



Table 1. Jeffreys guideline scale for evaluating the strength of evidence when two models are com- 
pared. 
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< 3 : 1 
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Inconclusive 


1.0-2.5 


~ 12 : 1 
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Significant 


2.5-5.0 


~ 150 : 1 


0.993 


Strong 


> 5.0 


> 150 : 1 


> 0.993 


Decisive 



Then, the quantity Bq^i provides an idea on how well model may fit the data when is com- 
pared to model 1. Jeffreys provided a suitable guideline scale on which we are able to make 
qualitative conclusions (see Table 1). In this work, we refer to positive (negative) values of 
i3o,i as the model being favoured (disfavoured) over model 1. 

The computation of the integral in Equation (2.3) is a very computationally demanding 
process, since it requires a multidimensional integration over the likelihood and prior. To 
carry out the exploration of the cosmological parameter space we first modify the CAMB 
code [23] in order to input an arbitrary shape of the primordial power spectrum^, then we 
incorporate into the COSMOMC software [24] a substantially improved and fully-parallelized 
version of the nested sampling algorithm called MultiNest [25, 26]. The MultiNest algo- 
rithm increases the sampling efficiency for calculating the evidence and allows one to obtain 
posterior samples even from distributions with multiple modes and/or pronounced degenera- 
cies between parameters. 



3 Primordial Power Spectrum Parameterisation 

3.1 Power spectrum reconstruction I 

First, we perform a reconstruction for the shape of the primordial spectrum. We parame- 
terise V-ji^k) with a specific number of bins, logarithmically spaced in k, and varying only 
each amplitude, denoted Ag^fc-. Throughout, we assume that most of the current relevant 
information is encompassed within the scales fcmin = 0.0001 Mpc~ and fcmax = 0.3 Mpc~ , 
where the combined WMAP-I-ACT data significantly improves the parameter constraints. 
Outside of these boundaries we take the spectrum to be constant with values equal to those 
at krain and fcmax respectively. We allow variations in the spectral amplitudes with a conser- 
vative prior As^ki £ [1, 50] x 10"^°. 

To maintain continuity between /c-nodes, a linear interpolation is performed such that the 



^A modified CAMB code version which allows a Node-based Parameterisation for the primordial power 
spectrum is a available at http://www.mrao.cain.ac.uk/software/ 



form of the power spectrum is described by 

Vnik) = I As,fc, k G {h} (3.1) 



^S.fcmin 
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k G {h] 


^S.fcmax 


n- ^ "'max 



and with hnear interpolation for fcmin < ki < k < ki+i < /cmax- 

We start our reconstruction by considering the base model is equivalent to the Harrison- 
Zel'dovich (HZ) spectrum {ug = 1) where the spectral amplitude VTi{k) = Ag is the only 
parameter, see Figure 1. The next model, (b), allows for two amplitudes located at k^i^ and 
^max to vary independently, thus emulating a tilted spectrum. We then add a third point (c) 
placed midway between the two existing nodes in (b) . This model mimics a degree of spec- 
tral running by allowing slight variations in the interpolated slopes between the three nodes. 
Since these amplitudes are independent of each other, however, there is no need to pick any 
particular pivot point as in the case of power-law parameterisation, hence providing more 
freedom in the shape of the spectrum. The presence of a turn-over in the global structure 
of the resultant spectrum, shown in Figure 1 (c), points to some deviation from a simple 
tilt. We might continue adding nodes in this fashion until some arbitrary accuracy of model 
fit is achieved, but always bearing in mind that the inclusion of new unnecessary nodes is 
penalised through the Bayesian evidence. We, then consider a fourth stage where the k-space 
is logarithmically split into three equally spaced regions, (d). At small scales the shape of 
the power spectrum is well constrained with tight error bars on each node, whereas on large 
scales the error bars tell us there is still room for new features (within the limited amount of 
information due to cosmic variance). Notice the increased error bars due to the addition of 
an arbitrary number of band-powers and correlations created between them. We also observe 
the evidence has dropped for four fc-nodes, therefore this stage seems to be a reasonable point 
to stop adding parameters in the reconstruction process. Figure 1 illustrates the correspond- 
ing form of the reconstructed spectra from the mean posterior estimates (with la error bars 
on its corresponding amplitudes), together with ID marginalised posterior distribution for 
the amplitude at each node and for each reconstruction. In each model we include the Bayes 
factor compared to the base model (HZ). 

The reconstructed spectra are assessed according to the Jeffreys guideline shown in Table 
1. The Bayesian evidence between the base model and the two-node model ;B2,i = +2.93ib0.30 
points out that the HZ is strongly disfavoured when compared to a tilted spectrum, in agree- 
ment with WMAP/ACT [2] results. The addition of complexity in the third stage provides 
more flexibility in the shape of the reconstructed spectrum. The evidence between the two- 
node and three-node model, i32,3 = +0.18 it 0.30, is too small to draw any decisive con- 
clusion, though the evidence marginally prefers the simple tilted spectrum. Although the 
reconstructed shape of the spectrum in the fourth stage is similar to the one obtained in the 
second stage, the four k-node model is penalised because of the inclusion of unnecessary infor- 
mation. Thus, the peak of the evidence at model (b) shows that the optimal reconstruction 
contains, surprisingly, just two nodes, as is shown in Figure 1. Therefore, according to this re- 
construction process, parameterisations such as the HZ and those containing more than three 
/c-nodes are hence disfavoured by current observations. At this point of the analysis, with 
fixed fe-node positions, our results are consistent with those obtained by [15], where accord- 
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Figure 1. Left: Reconstruction of the primordial spectrum modelled as piecewise linear between 
nodes with fixed wavenumber ki, along with mean amplitude values and their corresponding la error 
bars. On large scales the power spectrum is constrained by WMAP data, whereas at small scales 
ACT/LRG data yield tight constraints up to fc = 0.3 Mpc~^. Right: ID marginalised posterior 
distribution of the amplitudes Ai at each bin in each reconstruction. The top label in each panel 
denotes the associated Bayes factor respect to the base model (HZ) (a). 



ing to the Akaike information criterion, the preferred model is given by a two-node spectrum. 



To extract the global structure of the spectrum we have carried out a reconstruction 
process by placing nodes at particular positions in the A;-space. However, to localise features 



in A;-space, we may consider moving either back or forth the internal A;-nodes until we find 
their optimal position; we reconsider this option in an improved method in the next Section. 

3.2 Power spectrum reconstruction II 

In Section 3.1 we reconstructed the primordial spectruin using a standard binning process: 
fix k-node positions and vary only the amplitudes. We now consider a reconstruction of the 
spectrum where the internal k-node positions vary, as well as their amplitudes. 

3.2.1 Running-like spectrum 

In order to look for deviations from the simple power-law model, we consider a model with 
two fixed fc-nodes at sufficiently separated positions [A^min, ^max] > with varying amplitudes 
[^s.fcminJ ^s.fcniax]' ^^^ placc iusidc addltloual 'nodes' with the freedom to move around in 
both position ki and amplitude As^ki- Despite the simplicity of this approach, it covers a 
large variety of shapes for the primordial spectrum. The freedom of the position of the inter- 
nal /c-nodes allows us to localise the best position for a turn-over (if any) and the amplitudes 
are able to describe the global structure of the spectrum. 

Analogously to Section 3.1, we have maintained the same priors for the spectral am- 
plitude As^ki = [1)50] X 10^^'^, whereas on the fc-position we restrict to the physical prior 
log ki = [log A;min;log fcmax]- Hencc, for this type of nodal-reconstruction the spectrum is 
described by 



-^s,ki "^min ^ ki <. fcj+l < femax ("^-^j 

^s,fcniax '^ — kmax 



and with linear interpolation for fcmin < ki < A^max- 

The internal nodes generalise the spectral running by allowing slight variations in the inter- 
polated slopes between external nodes. Figure 2 illustrates the reconstruction of the shape 
of the primordial spectrum from the mean posterior estimates - with Icr error bars on the 
amplitudes - (left), along with the ID marginalised posterior distributions on the parame- 
ters used to describe the spectrum (right). On large scales, the reconstructed shape of the 
one-internal-node model (ki) resembles a similar spectrum to that obtained in Figure 1 (c), 
but now the probability distribution suggests a preferred turn-over position localised at the 
largest scales. A similar turn-over has also been identified using principle component analysis 
[17]. In the two and three-internal-nodes cases, it is interesting to note that at small scales 
the marginalised posterior peaks at scales where the combined WMAP/ACT constraints are 
improved; at these scales (0.1 < k [Mpc~^] < 0.14) CMB data now considerably overlap with 
measurements from SDSS DR7 LRG. This reduced power at small scales might be identified 
as a feature produced from a phase transition in the early universe [5]. Both, the two and 
three-internal models (middle and bottom panel of Figure 2), present a similar behaviour on 
the reconstructed spectra, also seen on the marginalised posterior distributions. 

3.2.2 Cut-ofF and Broken spectra 

For completeness, we consider the possible existence of a large-scale cut-off on the primordial 
spectrum. A possible motivation to consider this model has been discussed for instance by 
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Figure 2. Reconstruction of the primordial spectrum letting one internal fc-node (top) and two inter- 
nal fc-nodcs (bottom) move freely in both amplitude Ai and position ki. The right panel corresponds 
to the ID marginalised posterior distribution of the amplitudes and fc-node position in each recon- 
struction. The top label in each panel denotes the associated Bayes factor with respect to the base 
model (HZ) shown in Figure 1 (a). 



[27]. In order to perform the reconstruction for this particular case we fix an extremal node 
at fcmax with varying amplitude ^s,fcmax ^^^ l^t the cut-off scale kc vary across the prior 
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as well as its amplitude As,c- The form of the spectrum is described as follows: 



Vnik) = { As,c kc<k < fcinax 

^s.max "' r^ "-max 



(3.3) 



and with linear interpolation for kc < k < /cmax- 

Also, we consider a broken spectrum which might have been produced from phase transitions 
in the early universe. A similar broken spectrum, motivated by double or multiple field 
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Figure 3. Reconstruction of the large scale cut-off spectrum (top) and broken spectrum (bottom); 
their corresponding ID marginalised posterior distribution of the amplitudes Ai and node positions ki 
in each reconstruction (right). The top label in each panel denotes the associated Bayes factor with 
respect to the base model (HZ) shown in Figure 1 (a). 



inflation, has been considered by [5, 7]. This spectrum is obtained by placing two nodes fci 



and ^2 within [k 
such that 
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], and letting them move freely in amplitude ^s,fci and /c-position, 
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(3.4) 



and with linear interpolation for 



k-i < ko. 



The reconstruction of the cut-off and broken primordial spectra along with la limits of the 
marginalised distributions are shown in the left panels of Figure 3. Their corresponding pos- 
terior distribution in each parameter used to describe the spectra are illustrated in the right 
panels. The obtained best-fit parameters for the cut-off spectrum (top), show a preferred 
scale at which the power drops to zero with an upper limit log;^o ^c < —3.45 at 95% C.L. 
Our constraints on kc also show a significant likelihood at large scales, thereby disfavouring 
the presence of an abrupt cut-off. With respect to the broken model (bottom), on the other 
hand, the best fit parameters indeed predict a break in the primordial spectrum, located 
approximately at log^o ^ — —2.2. That could be an indicative of the existence of a phase 
transition, and it is similarly obtained in the two and three-internal-nodes models shown in 
Figure 2. 



In this section we have considered three types of spectra with different features: turn- 
over, large scale cut-off and broken spectrum (Figures 2-3). In each figure we have included 



the Bayes factor compared to the base model (HZ). According to the Jeffreys guidehne the 
one- internal- node spectrum, shown in the top panel of Figure 2, is significantly preferred over 
the cut-off and broken spectrum models, Bk^^k^ = +1-28 ± 0.30 and Bk^^kt = +1-88 it 0.30 
respectively. Even though the model with one-internal node is described by four parameters, 
when it is compared to the Harrison-Zel'dovich spectrum (with only one parameter), the 
Bayes factor B^^^i = -1-4.26 it 0.30 shows that the presence of a turn-over is strongly favoured 
by current cosmological information and significantly so when compared to the tilted spec- 
trum (with two parameters) i3fcj^2 = +1-33 it 0.30, see Figure 1. Therefore, the presence of a 
turn-over in VTi{k) plays an important role in explaining current observations. Notice that, 
in the bottom panel of Figure 2 the Bayesian evidence has dropped off, hence the reason we 
have stopped the addition of nodes in the reconstruction process. 

3.3 Power-law and running spectra 

We have considered, so far, a Vnik) shape reconstructed directly from data. For comparison 
we include the standard approach by assuming a power-law parameterisation in terms of a 
spectral amplitude Ag and a spectral index or tilt parameter Ug: 

/ I. \ "s-l 

Vn{k) = Ag i^-j , (3.5) 

where the spectral index is expected to be close to unity; /cq denotes the scale at which 
the tilted spectrum pivots, fixed to ko = 0.002 Mpc~ . We again assume the prior Ag = 
[1, 50] X 10^^'^ on the amplitude, together with the conservative prior Ug = [0.7, 1.2] on the 
spectral index. We find a mean value of ng = 0.963ib0.011 which confirms that our constraints 
are in good agreement with results from [1-3]. As a further extension we consider possible 
deviations from power-law by allowing the spectral index to vary as a function of scale ns{k). 
Then the primordial spectrum becomes 



rn{k) = Agl-] , (3.6) 

where Tirun is termed the running parameter and is expected to be nrun ~ for standard infia- 
tionary models. In order to minimise the correlation between ng and rirun we have considered 
a pivot scale of ko = 0.015 Mpc~^, as pointed out by [28]. We use the same priors as above 
on Ag and Ug, and the conservative prior rirun = [—0.3,0.3] on the running parameter. From 
the combined dataset we find the marginalised posteriors show a preference for a negative 
running parameter nrun = —0.026 it 0.015 and Ug = 0.968 it 0.011 for the spectral index, as 
expected by [1-3]. Figure 4 shows the marginalised posterior distributions for the parameters 
used to describe V-ji^k) and the obtained spectrum from mean posterior estimates of a simple 
tilted parameterisation with Ug (left) and including the running parameter nrun (right) re- 
spectively. In each panel we have included the Bayes factor compared to the base model (HZ). 

According to the Jeffreys guideline, present observations significantly prefer a simple 
tilted model when compared to a inodel which includes a running parameter by a factor 
Bns,nr:un = +1-19 ± 0.30. Similarly, a tilted spectrum is strongly preferred when compared 
to the HZ model: Bn^,! = -1-3.25 it 0.30. We also confirm the agreement between the simple 
tilted model and the two-fixed- noded spectrum through its Bayes factor, shown in Figure 
1 (b). An important point to emphasise is that the simple tilt and running model present 
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Figure 4. Reconstruction of tlic primordial spectrum assuming a simple tilted parameterisation 
with Us (top) and including the running parameter rirun (bottom); the coloured region denotes la 
error bands on the reconstruction. Right: marginalised ID and 2D probability posterior distributions 
for the power spectrum parameters; 2D constraints are plotted with Icr and 2cr confidence contours. 
The top label in each panel denotes the associated Bayes factor with respect to the base model (HZ) 
shown in Figure 1 (a). 



a significantly and strongly disfavoured Bayes factor, Bn^^^i = —1-01 ± 0.30, Bn^^^^.ki = 
—2.20 lb 0.30, compared to the reconstructed one-internal-node spectrum shown in Figure 
2. Thus a simple power-law parameterisation seems to be not enough to describe current 
cosmological observations, hence slight deviations of it should be taken into account. 

3.4 Modified Povi^er-law spectrum 

We have observed that models which present a turn-over at large scales are slightly preferred 
by the evidence. Based on this observation, we suggest the following phenomenological shape 
for the primordial power spectrum: 



VTi{k) = A, 






+ 1 



(3.7) 



In this particular parameterisation, assuming riv < 1, the parameter k^ determines the 
transition between a standard power- law model with red tilt {k ^ k^) to a blue tilt model 

{k<^ k^y. 



Vn{k) = A, 



riyj — l 






(3. 



11 



K) Sn.,1 = +4.65 ± 0.30 




0.0001 0.001 0.01 

-1, 





k[Mpc 



Figure 5. Reconstruction of the primordial spectrum from mean posterior estimates of the modified 
power-law parameterisation, along with la error bands (left). Right: marginalised ID and 2D proba- 
bility posterior distributions for the power spectrum parameters; 2D constraints are plotted with la 
and 2(7 confidence contours. The top label in each panel denotes the associated Bayes factor with 
respect to the base model (HZ) shown in Figure 1 (a). 

where the prior on n^ is similarly chosen to the spectral index Ug in the power-law pa- 
rameterisation: riy = [0.7,1.2]. We expect the constraints on the parameter k^ are mainly 
located on large scales, hence, for this extra-parameter we consider the following flat prior 
In /cv = [—15,-5]. The reconstruction of the shape of this spectrum along with the posterior 
distribution in each parameter are shown in Figure 5. We observe the constraints for the new 
tilt-parameter n^ = 0.955 it 0.014 are similar to those obtained from the power-law models, 
where the scale-invariant spectrum is ruled out at a high confidence level and the spectrum 
exhibits a red tilt at small scales. The marginalised posterior probability on k^ shows the 
existence of a blue-tilted spectrum on large scales In k^ < —8.1 at 95% C.L. Hence, the global 
shape for this spectrum presents a slight running behaviour with reduced power at both large 
and small scales, compared to the simple tilt parameterisation. 

The modified power-law parameterisation decisively rules out the HZ, since the Bayes 
factor between the models is Bn^,! = +4.65 ± 0.30. Moreover, even though the modified 
power-law model has an extra parameter compared to the simple power-law model, the data 
significantly prefer it with Bn^^ris = +1.40 ± 0.30. 

3.5 Lasenby & Doran spectrum 

The Lasenby and Doran model (LD, [6]) is based on the restriction of the total conformal time 
available in the entire history of a closed Universe. The primordial spectrum derived from 
this model naturally incorporates an exponential cut-off on large scales which might provide 
a possible explanation for the lower-than-expected CMB power spectrum at low multipoles. 

1/2 

On small scales, the relationship between V-j^ (k) and In k is linear, thus predicting a reduced 

1/2 

power at large k as compared to a simple tilted spectrum (for which InP^ (fc) versus In A: is 
linear). For further details about the LD model see, for instance [29, 30]. In this model, the 
resultant primordial spectrum from the inflationary phase depends upon just two parameters 
{bo,bi}. The parameter bo is mainly restricted by the number of e-folds A^ ~ 27rbf^ ~ 50, 
whereas 64 controls the initial curvature and must be negative and such that the dimension- 
less combination |64|/i~'^'^ is of order unity. Here fi is the mass of the scalar field, and this 
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Figure 6. Reconstruction of the primordial spectrum assuming a Lasenby & Doran model along 
with 1(7 error bands (left). Right: marginalised ID and 2D probability posterior distributions for the 
power spectrum parameters; 2D constraints are plotted with la and 2a confidence contours. The top 
label on the panel denotes the associated Bayes factor with respect to the base model (HZ) shown in 
Figure 1 (a). 



along with 6o and 64 control the magnitude of the inflaton field, and how long the inflationary 
period lasts. To compute the LD spectrum we refer our analysis to [19]. We have also chosen 
the priors based on the same paper: il^ = [—0.05, 10^^], b^ = [1,4], 64 = [—30, —1] x 10^^. 
Figure 6 shows the reconstructed shape of the primordial spectrum along with the posterior 
distribution in each additional parameter for this model; the constraints on the present bu- 
bble parameter are Ho = 69.4 it 1.4, whereas the number of e- folds is A^ = 50.6 it 4.3. From 
the top label of Figure 6, we observe the LD model is significantly preferred over the simple 
power-law parameterisation with a Bayes factor of ;BLD,ns = +1.69 ± 0.30 and decisive when 
compared to the HZ spectrcum: ;Bld,i = +4.94 ± 0.30. 

Indeed, the LD model has the largest evidence of all the models investigated, followed 
closely by the modified power-law spectrum. It should be noted, however, that the latter was 
constructed specifically to exhibit a turn-over on large scales, having already found that the 
data prefer such a feature. A fairer comparison would be between the LD model and the third 
most favoured model, namely the one-internal-node linear-interpolation model described in 
Section 3.2.1, since both of these models were proposed a priori. To this end, and as a check 
on our analysis, we use the best-fit LD model (shown in Figure 6) as the input spectrum to 
simulate an idealised CMB observation containing only cosmic-variance-limited noise. Fig- 
ure 7 (left panel) shows the resulting CMB temperature spectrum. We then analysed these 
simulated data using the one-internal-node linear interpolation model to reconstruct the pri- 
mordial power spectrum. Figure 7 (right panel) shows the resulting reconstruction (dotted 
line), which recovers well the shape of the input LD spectrum (solid line), except on the 
very largest scales, where there is little information in the simulated CMB data. Moreover, 
the reconstructed spectrum has a similar shape to the one obtained from real data using the 
one-internal-node model (see Figure 2). We may therefore understand the higher evidence 
for the LD model spectrum as resulting from its similar quality fit to the data, but requiring 
fewer free parameters than the one-internal-node linear-interpolation model. 

Finally, we note that, in the node-based reconstruction, the use of linear interpola- 
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Figure 7. Reconstruction of tlie Lasenby & Doran primordial spectrum based on the binning format 
with linear interpolation described in Section 3.2.1. We have assumed an idealised CMB spectrum 
with limitation only due to cosmic variance (left). Right panel shows the reconstructed spectrum in 
the binning format together with the LD input spectrum (solid line). 

tion between the nodes may seem crude. It is straightforward to generalise the node-based 
approach to more sophisticated interpolation schemes, but this may not always yield better 
results. In Appendix A, we illustrate this point by reanalysing the simulated CMB data using 
a cubic spline interpolation through the nodes, thus allowing one to reconstruct a smooth 
shape for the primordial spectrum, but one that is less satisfactory than that obtained using 
linear interpolation. 

4 Discussion and Conclusions 

In this paper we have attempted to fit an optimal degree of structure for the primordial 
power spectrum of curvature perturbations using Bayesian model selection as our discrimi- 
nating criterion. We have modeled the spectrum as a linear interpolation between a set of 
'nodes' with varying amplitude and /c-position. We have also explored different parameteri- 
sations of the primordial spectrum which include: a power-law parameterisation with both 
tilt and running parameter, a modified power-law spectrum and the Lasenby & Doran (LD) 
model. 

All the considered models have in common the standard cosmological parameters: r^b^^, 
f^DM^^, ^, ''", as well as the secondary parameters: Asz, Ap,Ac. Thus, priors on these param- 
eters remained the same in each model. The best-fit values for these standard parameters 
are consistent with those obtained using the concordance 6-parameter model with power-law 
primordial spectrum. We show, in Figure 8, ID marginalised posterior distributions for the 
cosmological parameters of each of the preferred models. We observe that the values of the 
standard parameters remain well constrained despite the addition of extra freedom on the 
shape of the primordial spectrum, although the constraints resulting from the HZ spectrum 
clearly differ from the others. Note also that the constraints on the parameters corresponding 
to the LD model are slightly tighter than the rest of the models. 

We have considered wide-enough priors in our analysis, such that they do not interfere 
with the inferred parameter values. We used priors on the amplitudes of Ai = [1, 50] x 10^^" 
and on spectral indices of ris = n^ = [0.7, 1.2], while the parameters describing the /c-space 
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Figure 8. ID marginalised posterior distributions for the primary (top) and secondary (bottom) 
cosmological parameters, for each corresponding model listed on the bottom-right box. 



have physical priors restricted by [A;mini ^max]- We also compute the Bayesian evidence for a 
wider prior range of Ug = n^ = [0.5, 1.5] and 64 = [—50, — 1] x 10^^, to illustrate the robust- 
ness of a model over small variations of the prior range: 



B„,,,i = +2.25 ± 0.30 (wide priors) 
Bn^^i = +4.24 ± 0.30 (wide priors) 

Sld,i = +4.47 ± 0.30 (wide priors) 

We observe that even when wider priors are considered, the HZ model is strongly disfavoured 
when compared to n^ and LD models. Similarly, the simple tilted model is still significantly 
disfavoured. 

To summarise the analysis, in Figure 9 we plot the reconstructed spectra for the pre- 
ferred selected models together with their corresponding Bayesian evidence. It shows that 
the HZ spectrum is decisively excluded as a viable model to describe Vii{k). The preferred 
model given current observations is provided by the LD model followed by a modified power- 
law version. We have found that the power-law parameterisation, including either cases Ug, 
and ng + rirurn are both significantly disfavoured. The presence of a turn-over at large scales ^ 
and the reduced power at small scales seem to provide an important contribution on choosing 
the best-fit model through its Bayesian evidence. 



At the largest scales, the addition of tensor perturbations might considerably affect the shape of the 
primordial spectrum. Hence, we consider this possibility in a more detailed future work. 
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Figure 9. Comparison of the primordial power spectra for the preferred models along with their 
Bayesian evidence. We also include the maximum likelihood £,nax for a model with number of 
parameters Npar- Each Bayes factor is compared respect to the one-node model (HZ). 
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A Cubic spline interpolated spectrum 

In this appendix we illustrate how the type of interpolation in our node-based approach can 
influence the reconstruction of the primordial spectrum. We use the same example shown 
in Figure 7 but now we use a cubic spline to interpolate through the A:-nodes. From Fig- 
ure 10 we note that the spectrum used as an input lies well outside the error bar on the 
reconstruction at low /c-values. Therefore, the spline fails to recover the input spectrum cor- 
rectly, contrary to the linear interpolation where the recovered spectrum, shown in Figure 
7, is more representative, with the true spectrum lying comfortably within the error-bars on 
the reconstruction at all /c- values. This is mainly because a function with rapidly changing 
higher derivatives, such as the input primordial spectrum used here, is less accurately ap- 
proximated by higher order polynomials. In particular, the requirement of continuous first 
and second derivatives, combined with the tight constraints at small and intermediate length 
scales, leads to a significant overestimation of the power at the less well constrained region 
at the largest scales. Hence, in this case, the linear interpolation describes the shape of the 
primordial spectrum better than the spline. 

For comparison with the results presented earlier, we also use the cubic spline to perform 
similar node-based reconstructions of the spectrum to those shown in Figures 1 and 2. The 
low number of bins used to describe the global structure of the spectrum yield to similar 
shapes by using both interpolation methods. The reconstructed spectra for three and four 
bins, along with one and two internal A:-nodes, are plotted in Figure 11 using the spline 
interpolation. 
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Figure 10. Reconstruction of the Lasenby & Doran primordial spectrum based on the binning 
format with cubic spline described in Section 3.2.1. We have assumed an idealised CMB spectrum 
with limitation only due cosmic variance (left). Right panel shows the reconstructed spectrum in the 
binning format together with the LD input spectrum. 
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Figure 11. Reconstruction of the primordial spectrum using the cubic spline. Top panel resembles 
plots (b) and (c) shown in Figure 1, whereas bottom panel the reconstruction for the models shown 
in Figure 2. To the right of each reconstruction we plot the ID marginalised posterior distribution of 
the amplitudes Ai and fc-node position ki. The top label in each panel denotes the associated Bayes 
factor with respect to the base model (HZ) shown in Figure 1 (a). 
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